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Abstract 

Working within an exactly solvable 3 level model, we discuss an extension 

of the Random Phase Approximation (RPA) based on a boson formalism. 

A boson Hamiltonian is defined via a mapping procedure and its expansion 

truncated at four-boson terms. RPA-type equations are then constructed and 

solved iteratively. The new solutions gain in stability with respect to the 

RPA ones. We perform diagonalizations of the boson Hamiltonian in spaces 

containing up to four-phonon components. Approximate spectra exhibit an 

improved quality with increasing the size of these multiphonon spaces. Special 

attention is addressed to the problem of the anharmonicity of the spectrum. 
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I. INTRODUCTION 



The most commonly used microscopic approach for the study of collective vibrational 
states in many-fermion systems is the Random Phase Approximation (RPA) In this 
theory the lowest collective excitations result from the action of phonon operators Q\, on a 
state \RPA) which is defined by the condition that Q U \RPA) = 0. This state represents the 
ground state of the system. It is a distinctive feature of RPA that multiphonon states, i.e. 
states obtained by repeated actions of phonon operators on the ground state, are eigenstates 
of the Hamiltonian with energies forming a harmonic spectrum. The existence of states 
which can be approximately described as corresponding to the multiple excitation of low- 
lying and/or high-lying phonons is well established in atomic nuclei. However, deviations 
from the harmonic picture are also observed and their influence on several processes has 
been analysed 0. 

In a standard derivation of the RPA equations, a crucial point is represented by the 
so-called Quasi-Boson Approximation (QBA). This is a rather crude approximation which 
causes the Ql operators to behave as boson operators in spite of their (composite) fermionic 
structure. Overcoming this approximation has represented the starting point of many at- 
tempts aiming at improving RPA p|-p2|. One of the line of research in such a context 
has been based on a reformulation of the whole theory in a boson formalism [|T^-^]. In 
other words, the Q' u operators have been defined from the beginning in terms of true boson 
operators and all the fermion operators of interest have been replaced, via a mapping proce- 
dure, by their boson images. The RPA-type equations that one constructs in this formalism 
depend, of course, on the degree of expansion of the boson Hamiltonian. Truncating this 
expansion at the lowest order, i.e. at two-boson terms only, gives the boson counterpart of 
RPA. Including higher-order terms in the boson image of the Hamiltonian provides a natural 
way to reach a higher level of approximation. Besides that, the inclusion of these terms has 
another important effect: it leads to a coupling among multiphonon states. States which 
result from a diagonalization in a m-phonon space are therefore superpositions of zero-, 
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one-,... , m-phonon states. Such a diagonalization is expected to lead to a further improved 
degree of approximation as well as to cause anharmonicities in the spectrum. 

Calculations in this boson formalism have been performed in the recent past for atomic 



nuclei |19|]2(| considering a Hamiltonian truncated at four-boson terms and diagonalizing it 
in the space of one- and two-phonon states. The resulting anharmonicities have not been 
found large, especially in 208 Pb. In particular, the anharmonicity associated with states 
whose main component is a Double Giant Resonance has been found of the order of a few 
hundred KeV. This is certainly related to the fact that RPA gives a good description of 
Giant Resonances, especially in heavy closed shell nuclei. 

In metallic clusters, a strongly collective state is known to exist, the dipole plasmon, 
which corresponds to the oscillation of the delocalized electrons of the cluster against the 
positively charged ions. The experimental evidence for states corresponding to the double 



excitation of the plasmon has not been confirmed p3[ . From the theoretical point of view 
the situation is also quite unclear. From one hand, in a purely harmonic spectrum for 
the multiple excitation of the plasmon has been predicted. On the other hand, by using the 



same approach as in |fT5| , p(l| huge anharmonicities in the two-plasmon states have been found 



[25|. An important difference with respect to the case of atomic nuclei is that in metallic 



clusters the two-body interaction is very long ranged. This is probably the main reason why 
the RPA ground state is very different from the Hartree-Fock one and the RPA backward 
amplitudes are quite large. Of course, this may cause that the same level of truncation in 
the boson expansion is not adequate both in the case of nuclei and of metal clusters. 

In principle, configuration mixing calculations can give a clear, model independent, in- 
dication on the existence of such two-phonon states and on their degree of anharmonicity. 
Unfortunately, since the states one looks for are quite high in energy, the number of config- 



urations required to get stable results is huge. In |26[] such a study has been performed for a 
very simple case: two interacting electrons moving in a uniform positive charge distribution. 
This is a kind of precursor of a metal cluster in the jellium approximation and allows for 
a numerically exact calculation. Important deviations from the harmonic limit have been 



found. More specifically, in addition to an almost perfectly harmonic vibrational band based 
on the ground state, the other states can be classified in other bands with a much less degree 
of harmonicity. 

Aim of this paper is to shed some light on the limits of the approach adopted in |2D| and 
by applying it to a 3-level solvable model E7\. The analysis we are going to present 



is very similar to that in [^TJ where a 2-level model was considered and the parameters 
were adjusted in such a way to mimick the multiple excitation of a Giant Resonance. The 
3-level model is, of course, richer. In particular, since there are two single particle states 
above the Fermi surface (particle states) and one below (hole), two different elementary p-h 
configurations and, correspondingly, two different phonons can be excited. Therefore, one 
can better simulate the situation encountered in nuclei which generally present one high- 
lying and one low-lying collective modes for each multipolarity. Also, matrix elements of 
the interaction connecting a particle-hole state with a two-particle one can be included in 
a natural way. These terms are present in a generic two-body interaction and are very 
important since they couple states having numbers of phonons differing by one. 

The paper is organized as follows. In Sec. II, we will describe the model and the 
formalism. In Sec. Ill, we will analyse the results. Finally, in Sec. IV, we will review the 
content of the paper and draw some conclusions. 



II. THE MODEL AND THE FORMALISM 

The model p7| consists of three 2f2-fold degenerate single-particle shells which are oc- 
cupied by 2Q particles. In the absence of interaction, then, the lowest level is completely 
filled while the others are empty. This state, the "Hartree-Fock" (HF) state of the system, 
is denoted by |0). A single-particle state is specified by a set of quantum numbers (j,m), 
where j stands for the shell (j =0,1,2) and m specifies the 2Q substates within the shell. The 
creation and annihilation operators of a fermion in a state (J, m) are defined by aj m and 
djm, respectively. 
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Let us consider the operators 

2Q 



K ij = E a im a jm = 0, 1, 2). (2.1) 
m=l 

These operators satisfy the Lie algebra of the group SU(3) 

[K{j, K ki ] = S jk Ku - 5uK h j. (2.2) 

The Hamiltonian of the model is written in terms of the generators only and contains 
up to two-body interactions. Its form is 

H F = J2e(i)K u + V x (i,j)K i0 K oj 

i=l,2 »J=1,2 

+\ E V v (i,j)(K io K jQ + K oj K oi ) 

+ V y {i,j,k)(K iQ Kj k + K kj K oi ), (2.3) 

i,j,k=l,2 

with real coefficients. The eigenstates of Hp can be constructed either by using the properties 
of the SU(3) algebra or by diagonalizing it in the space 

nina) = 7 =^(K 10 )-(^ 20 ) n2 |0) 1 , (2.4) 

V 7 ^ 1 " 2 > 0<n 1+ n 2 <2n 

where M ni n 2 are normalization factors. 

In order to introduce the boson formalism, let's define the space 

B = {\n x n 2 ) = -J=(&t r(& t no) ] ; (2 _ 5) 



\Jn x \n 2 \ 



0<ni+n 2 <2n 



where the operators b\ obey the standard boson commutation relations 

M]]=%, [bi,bj) = (2.6) 

and |0) is the vacuum of the b^s operators. A one-to-one correspondence exists between the 
states of F and B, the boson operators b\ playing the role of the excitation operators K i0 
and the boson vacuum |0) replacing the HF state |0). 

The mapping procedure to construct boson images of fermion operators is the same dis- 
cussed in previous works (see, for instance, Ref. [p8fl ) and, due to the orthonormality of both 



sets of states \ni, n 2 ) and \rii, 77.2), it is simply based on the requirement that corresponding 
matrix elements in F and B be equal. The procedure is, therefore, of Marumori-type. We 
refer to Ref. ]2S[ for more details. Here, we simply say that, in correspondence with the 
Hamiltonian Hp ( |2.3| ), we introduce a hermitian boson Hamiltonian Hb which, in the most 
extended version, has the form 



H B = a + X PiQi + h.c.) + X nM b j + X M^b] + h.c.) 

i ij i<j 

+ X X e m(blb]b k + h.c.) + X X $ijkib\b]b k bi. 

i<j k i<j k<l 



(2.7) 



with the coefficients depending on the parameters e(i), V x ,V v ,V y of Eq.(^l| 
Let us now introduce the operators 

Ql = X*^!-X>fV 

i i 

and let the state |\l/o) satisfy the condition 

Q v |*o) = 0. 



(2.8) 



(2.9) 



By using the equations of motion method |29[ one finds that the amplitudes X and Y are 
solutions of 



A B 

-B* - A* 









f ^ 








J 


V(«") 

V J 







(2.10) 



where 



= (*o|[6 ij [^B,6j]]|*o), 

S fi = -(*o|[6i,[^B,6i]]|*o) 
and u;^ are the energies of the excited states 

M = Qj|*o). 



(2.11) 
(2.12) 



(2-13) 



As anticipated in the Introduction, the form of Eqs. fl2.1C|) is strictly related to the degree of 
truncation of the boson Hamiltonian. In the hypothesis that Hb contains up to two-boson 



terms, the double commutators in Eqs. (|2.11 ) and (|2.12|) are just numbers which, therefore, 



are also the values of the matrices A and B. This is the simplest case which can be realized 
in this formalism and represents the boson counterpart of the standard RPA. This degree of 
approximation can be improved by introducing a Hamiltonian with higher-order terms like, 
for instance, ( |2.7[ ). In this case the double commutators are operators. In order to calculate 
their expectation values in |\l/o) > as required in Eqs. (|2.11| ) and ( p,12|) , one can express the b 



and operators in terms of Q and by reversing Eq. (|2.8| ) (and its conjugate) and use the 
fact that the ground state |\l/o) is defined as the vacuum of the Q's. This procedure gives, 
however, matrices A and B which depend on the X and Y amplitudes and, consequently, 
equations of motion ( |2.10|) which are non linear. In what follows this non linear extension 
of RPA will be called ERPA. 

Having determined the X and Y amplitudes within RPA or ERPA, one can express the 
Hamiltonian Hb in terms of the operators Q and . In the case of RPA, namely when the 
boson Hamiltonian (|2.7|) is truncated at two-boson terms only, Hb can be rewritten simply 



as 



Hb = E + J2^ ] QIQu, (2.14) 



where are the energies solutions of the RPA equations Q2.10I ). This Hamiltonian obvi- 



ously doesn't mix states with different phonon numbers and so its eigenstates are pure zero-, 
one-, m-phonon states. For a higher-level truncation in the boson Hamiltonian, like for 
instance that of Eq. (|2.7| ), Hb acquires instead the more general form 



H B = E + H 10 (Q^ + h.c.) + H U Q^Q + H 20 (Q^Q^ + h.c.) + 
H 2 i{Q ] Q ] Q + h.c.) + H^(Q ] Q ] Q ] + h.c.) + 

H 22 Q ] Q ] QQ + H^{Q ] Q ] Q ] Q + h.c.) + H m {Q ] Q ] Q ] Q ] + h.c), (2.15) 

(for simplicity, we have dropped all the indices) where the Hij coefficients are functions of 
X and Y. Also in this case, as in RPA, the term H 2 o as well as the non-diagonal terms Hn 
vanish, as can be easily shown using the fact that the X and Y amplitudes are solutions of the 



ERPA equations ( |2.10|) . The remaining terms of ( |2.15 ), however, mix states with different 



phonon numbers so that the eigenstates of the full Hamiltonian become combinations of these 
states. This fact introduces an evident difference with RPA since the energies which result 
from the ERPA equations are not eigenvalues of the boson Hamiltonian in the phonon space 
as it is in the case of RPA (where they provide the excitation energies of the one-phonon 
eigenstates). In the next section, we will show the results obtained by diagonalizing Q2.15|) 



in different bases, containing up to two-, three- and four-phonon states, and compare them 
with the RPA and the exact ones. 

III. RESULTS AND DISCUSSION 



Let us first of all fix the parameters entering in the Hamiltonian ( |2.3| ). The energies 
e(l) and e(2) have been chosen equal to e and 2.5e, respectively. In order to simplify the 
analysis, the coefficients V x (i,j), V v (i,j), V y (i,j,k) have been assumed independent of the 
indices i,j, k and proportional to one parameter, \. More precisely, we have fixed V x = — x, 
Vv — \Xi Vy — ~ §X as i n Ref- 0- Both e and x are parameter expressed in units of energy. 
All calculations have been performed for a system of 2Q = 10 particles. 

We begin our analysis by examining the limits of validity of the truncation in the boson 
Hamiltonian (|2.7|). To do that we compare the exact spectrum of the fermion Hamiltonian 



1|) with that obtained by diagonalizing Hb (|2.7| ) in the full boson space. The results are 
shown in Fig. 1 where we report, as functions of the strength r = the energies of those 
states which, at zero strength, are pure lp-lh or 2p-2h states. For later use we also label 
these states with the symbols \u) and \v1U2) meaning by that the states which, in the limit 
r — > 0, are the one-phonon = Ql\ty )) and two-phonon (\u1U2) = Q^QlJ^fo)) states of 
the RPA formalism. The eigenvalues of the boson hamiltonian (dashed lines) are found in 
good agreement with the exact ones in the whole range of variation of r. In the following 
we will take the spectrum of the truncated hamiltonian as the reference one. 

In Fig. 2, we show the energies resulting from the RPA and ERPA equations together 
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with the reference spectrum. For small values of the strength, RPA and ERPA give almost 
identical results, in good agreement with the reference ones. Differences become really visible 
around r = 0.28, where RPA undergoes a collapse. In this region ERPA solutions remain 
instead stable and close to the reference ones. On the basis of the comments at the end of 
Sec. II, however, we notice that these ERPA energies can be considered excitation energies 



of one-phonon states only in the hypothesis that all the terms of ( p.!5|) but E$ and Hn play 
a minor role. 

Although not clearly visible at a first glance, the spectrum of Fig. 1 shows marked 
anharmonicities. In order to quantify these, in Fig. 3 we plot the ratios 

where by E v and E VlVa we mean the exact energies of the states \v) and \v1V2) defined in 
Fig. 1. One notices a well different behaviour of the three ratios plotted. In the case of the 
state 1 11) the ratio is seen to grow rapidly for increasing r. In the case of the state |22), 
R V1V2 remains instead almost exactly zero therefore showing that the corresponding exact 
state can be fairly well described as a pure two-phonon state. The remaining "mixed" state 
1 12) exhibits only a moderate anharmonicity. These results resembles that of Ref. [2B| where 



a realistic two-electron system was examined. Indeed, in both cases, the spectrum is found 
to exhibit some levels with a clear harmonic nature mixed, however, to other levels which 
do not display at all this nature. This mix of states makes of course complicate to provide a 
clear cut answer to the question of how much harmonic is the spectrum of these many-body 
systems, the answer strongly depending on which levels one is looking at. 

The existence of the anharmonicites evidenced in Fig. 3 represents an evident limit to 
the harmonic picture of RPA and points to the need of diagonalizing the full Hamiltonian 
( |2.15| ) in a multiphonon space in order to reproduce the exact spectrum properly. We have 
performed three separate calculations by diagonalizing Hb in spaces including up to two-, 
three- and four-phonon configurations. The results are shown in Fig. 4 together with the 
reference ones. The internal structure of the collective phonons entering in the basis states 

9 



is the one obtained in ERPA, i.e. the X and Y amplitudes are solutions of Eqs. ( [2.10 ) where 



the full Hamiltonian (|2.7| ) has been used. Keeping for simplicity the same labelling of states 
used in Fig. 1, we can say that, as expected from the previous discussion, for the states |1), 
1 2) and 1 22) the results are almost independent of the basis and in very good agreement with 
the exact ones. For the state |11) we find instead a marked dependence on the basis, the 
results approaching more and more the reference ones when the space is enlarged. Going 
from the smaller space to that including up to three phonons the agreement is quite good for 
r < 0.25. For a higher strength the space needs to be further enlarged up to four phonons. 
A somewhat similar behaviour, although less pronounced, is seen for the state |12). 

As a final result, in Fig. 5 we show the ground state energies calculated within RPA 
and as they result by diagonalizing the Hamiltonian ( |2.15|) in multiphonon spaces. All these 
results are compared with the exact ones and with those obtained by diagonalizing the 
boson Hamiltonian ( |2.7|) . It is interesting to see that the good agreement between fermion 
and boson results already found in Fig. 1 for the excitation energies is confirmed also in 
this case. As expected, RPA overestimates the exact values. All the diagonalizations in 
multiphonon spaces lead instead to an underestimation of these values showing, however, a 
rapid convergence to the exact values with increasing the size of the spaces. 



IV. CONCLUSIONS 

In this paper we have analysed some of the lowest excited states of the spectrum of 
a solvable 3-level model, namely those that in RPA would be described as one and two- 
phonon states. We have worked in a boson formalism. As a preliminary step, then, we have 
constructed a boson image of the fermion Hamiltonian whose expansion has been truncated 
at four-boson terms. The procedure followed in such a derivation has been of Marumori- 
type. The quality of the boson Hamiltonian has been tested by comparing its eigenvalues 
with the exact ones. Within the considered range of variation of the interaction strength 
the agreement between fermion and boson energies has always been found quite good. 
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By making use of this boson Hamiltonian we have constructed RPA-type equations and 
solved them iteratively. The new solutions have gained in stability with respect to the RPA 
ones and, in particular, around the RPA collapse point the new energies have exhibited a 
good agreement with the exact ones. This extension of RPA, while introducing corrections to 
the Pauli Principle violations present in RPA, naturally leads to a Hamiltonian which mixes 
states with a different number of phonons. We have performed diagonalizations in spaces 
containing up to two-, three- and four-phonon states and observed an improved quality of 
the approximate spectra with increasing the size of the spaces. 

Special attention has also been addressed to the problem of the anharmonicity of the 
spectrum. This has been found relevant for the state which, in RPA, is described as the 
double excitation of the lowest one-phonon state. On the contrary, for the other states, a 



less pronounced anharmonicity has been found. These findings agree with those of Ref. [26 



and point to the necessity of considering together all the possible elementary excitations of 
a many-body system when discussing the anharmonicity of its spectrum. 
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FIGURES 

FIG. 1. Excitation energies of the states \u > and \ v\V2 > which, at zero interaction strength, 
are pure lp-lh and 2p-2h states respectively. The exact energies (full lines) are compared with the 
energies obtained diagonalizing the hamiltonian Hb in the full boson space (dashed lines). The 
energies are shown as functions of the parameter r. 

FIG. 2. Excitation energies, as functions of r, calculated within RPA (dotted lines) and ERPA 
(dashed lines) compared with the corresponding reference energies (full lines). 

FIG. 3. The ratios R Ul u 2 (Eq.fljP])) calculated for the states |11 >, \12 > and |22 > as functions 
of r. 

FIG. 4. Excitation energies, as functions of r, obtained diagonalizing Hb in spaces including 
up to two- (dotted lines), three- (dot-dashed lines) and four- (dashed lines) phonon configurations 
compared with the reference excitation energies (full lines). 

FIG. 5. Ground states energies, as functions of r, obtained in the exact fermionic calculation 
(full line), in the reference calculation (dashed line), in RPA (+ symbols), in the diagonalizations up 
to two phonons (left triangles), up to three phonons (diamonds) and up to four phonons (circles). 
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